clc
tic
for jj=1:5
likelihood(data,setup.initial_parameter,setup)
end
toc